Co-operative Kondo Effect in the two-channel Kondo Lattice 
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We examine how the properties of a single-channel Kondo lattice model are modified by additional 
screening channels. Contrary to current wisdom, additional screening channels appear to consti- 
tute a relevant perturbation which destabilizes the Fermi liquid. This instability occurs involves two 
stages. When a heavy Fermi surface develops, it generates zero modes for Kondo singlets to fluctuate 
between screening channels of different symmetry, producing a divergent composite pair susceptibil- 
ity. Additional screening channels couple to these divergent fluctuations, promoting an instability 
into a superconducting state with long-range composite order. We discuss possible implications for 
heavy fermion superconductivity. 
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I. INTRODUCTION 

One of the remarkable properties of localized magnetic 
moments is their ability to transform the electronic prop- 
erties of their host. rTihese effects are dramatic in heavy 
fermion compoundsB'EI Since the mid-seventies, several 
hundred heavy fermion compounds have been discovered, 
characterized by a dense lattice of magnetic rare-earth, 
or actinide ions immersed in a conducting host. These 
materials by-pass the normal development of ordered an- 
tiferromagnetism to form a new kind of electron fluid. 
The resulting metallic state contains quasiparticles with 
effective masses up to a thousand times greater than a 
bare electron. For example, in CeCu$l the presence of 
only 14% Cerium in the copper host increases the effec- 
tive mass of the electrons by a factor of 1600. 

In a small handful of heavy fermion compounds, the 
heavy electron fluid becomes superconducting^ Local 
moments, usually harmful to superconductivity actually 
participate in this superconducting condensation pro- 
cess and a significant fraction of the local moment en- 
tropy is quenched as part of the condensation process. 
In UBeis for example, the spin-condensation entropy is 
about 0.2fcs?n2 per spin.cl One of the great challenges is 
to understand how microscopic order parameter in these 
systems involves the spin operators of the local moments. 

The concept of "composite pairing", where a Cooper 
pair and local moment form a bound-state combination 
that collectively ccoidejises may provide a way to ad- 
dress this problem.trEiJ A composite "triplet" involves 
a bound-state between a spin and singlet Cooper pair, 

A t (x) = (^ N - 2 \yj l (x)^(x)S(x)\^ N ), (1) 

but a composite singlet involves a triplet and a spin-flip, 

A s (x) = (y N _ 2 \vj n (x)^ 2l (x)S + (x)\y N ). (2) 



where "1" and "2" refer to two conduction electron chan- 
nels. Such composite order parameters were originally 
considered in the context of odd- frequency pairingpQ but 
more recent work has emphasized that composite order 
may co-exist with BCS pairing in cases where-.the spin 
plays a central role in the condensation process Ba Unfor- 
tunately, we know very little about how such composite 
pairing might come to pass. A divergent composite sin- 
glet susceptibility is known to occup-.in the symmetric 
two-channel Kondo impurity model,El and more recent 
studies suggest that a large composite susceptibility may 
persist into the two-channel Kondo latticeO 

In this paper we introduce a model for heavy fermion 
behavior where the local moments couple to a single con- 
duction band via two orthogonal scattering channels. We 
find that when two scattering channels of the same par- 
ity share a common Fermi surface, constructive interfer- 
ence develops between the channels. The scattering of 
electrons in the Kondo effect is described by an SU(2) 
matrix Vr, (r = 1,2) associated with each channel. A 
key result of our paper relates the composite order to 
the gauge invariant interference term between these two 
matrices: 



V f 2 Vi = - 



F A 

-At Ft 



where 



F = V>*lO"02 • S, 

A = ipi(ia y )<np2 ■ S, 



(3) 



(4) 



represent the singlet composite order in the particle-hole, 
and particle-particle channels respectively, and J\ and 
J 2 describe the Kondo coupling constants in the two 
channels. In a single impurity, the Kondo effect in the 
stronger ehaanel, suppresses any Kondo effect in weaker 
channels .OEJ A key feature of our lattice mechanism, 
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FIG. 1. Mean-field phase diagram for a two-dimensional 
two-channel Kondo lattice with an s-wave and d-wave inter- 
action channel. n c denotes the filling of the conduction band. 
J2/J1 is the ratio of coupling constants in the two channels. 
The phase diagram was computed for a tight-binding model, 
keeping max(Ji, J 2) = At, the band-width. "K.I " denotes 
"Kondo insulator" phases, that exist at half-filling, but which 
undergo a superconductor-insulator transition at a critical 
value of J2/J1. The lightly shaded region is dominated by 
composite pairing, and there is a gap for quasiparticle exci- 
tations. Along the dotted lines, the conduction electrons are 
entirely unpaired, and a pure composite pair condensate is 
formed. In the darkly shaded regime, Cooper and compos- 
ite pairs co-exist, and a gap less anisotropic superconductor is 
formed. 



is that channel interference co-operatively enhances the 
Kondo effect in the weaker channel, driving the develop- 
ment of composite pairing for arbitrarily weak second- 
channel coupling. 

The development of phase coherence between the two 
channels is signalled by the condensation of composite 
singlet pairs at a new temperature scale 



(5) 



where Tri and Tk2 are the Kondo temperatures of the 
primary and secondary channels respectively. The un- 
derlying gap symmetry of the quasiparticles in this new 
superconducting phase reflects the interference phenom- 
ena and is is given by a product of the form factors $ik 
and <£>2k fr° m each channel: 



A k = A $i k $2k- 



(6) 



In the typical composite paired state, composite pairs 
co-exist with Cooper pairs, as envisaged HVithe works of 
Bonca and Balatsky and also Poilblanc.Q'cl One of the 
novel features of this mechanism, is that it permits both 
gapless, and gapfull anisotropic superconductivity. In the 
region where the coupling constants J\ and J2, for the 



two channels, are of comparable strength, the nodes in 
the excitation spectrum gravitate to the center of the unit 
cell, where they mutually annihilate to produce a gapped 
phase (Fig. [j]., Fig. [?].). In this phase, the BCS order 
parameter is very small, and actually vanishes along a 
line in the phase diagram. 

At half filling, .a_Kondo lattice generally forms a Kondo 
insulating phaseo With the addition of a second-channel 
coupling, we find that at a critical ratio of coupling con- 
stants, there is a second-order transition from the Kondo 
insulating phase into a pure condensate of superconduct- 
ing pairs. This leads to a phase diagram, where a first 
order line representing the Kondo insulator terminates at 
a superconducting-insulating transition, as illustrated in 
Fig. d). 

An brief description of our work in this area has already 
been publishedJIj This paper is intended to provide a 
detailed account and discussion of the co-operative two- 
channel Kondo effect. 



II. MULTI-CHANNEL SCATTERING EFFECTS 
IN INTERACTING KONDO LATTICES 

The classical approach to heavy fermion physics in- 
volves local moments which couple exclusively to con- 
duction electron states with the same local f-symmetry. 
This assumption derives from the observation that spin- 
exchange between the conduction electrons and the local 
moments occurs predominantly via hybridization in the 
f-channel. ■— ■ ,— ■ 

However, more careful considerations suggestliZrtL£l 
that electron-electron interactions can cause new spin- 
exchange channels to open up between a local moment 
and the conduction sea. There are several mechanisms 
by which these new spin exchange channels can open up, 
including 

• Vicinity to a quantum critical point .0 



• Interactions in the conduction sea.ESE; 

• Intra-atomic Hunds interactions^ 

The first mechanism, identifiedJpng ago in a little-known 
paper by Larkin and Melnikovll3 may be particularly im- 
portant for heavy fermion systems which lie at the brink 
of magnetism. Larkin and Melnikov studied the single 
impurity Kondo effect in the vicinity of a magnetic quan- 
tum phase transition, where the local local moment po- 
larizes electrons at increasingly greater distances. The 
critical magnetic order thereby induces the spin to scat- 
ter electrons in a large number of angular momentum 
channels up to a maximum value l ~ kp^, where £ is 
the spin correlation length. The large screening cloud 
causes the matrix element for spin exchange to become 



J -> Jk,k' = ^x(k - k'), 



(7) 
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FIG. 2. Magnetic moment in an interacting environment. 
Localized electron at center of plaquet hybridizes in the 
d^-channel with nearby atoms. The on-site interaction at 
each atomic site U c is taken to be far larger than the electron 
band-width t. 

where x is the strongly momentum dependent suscep- 
tibility of the magnetic host. When decomposed into 
partial wave states, they found that this led to a Kondo 
coupling to electrons in all channels with angular mo- 
mentum I < l — kp^. 

More recent work has made it clear that new spin ex- 
change channels open up whenever charge fluctuations, 
are suppressed by interactions in the conduction sea£3't£l 
Consider the situation shown in Fig g. , where a local mo- 
ment hybridizes with nearby orbitals in a a d-channel. 
The spin-exchange between the local moment is written 

Hi = J(S • ^ da a a ^ d p), (8) 

where 

*L = ^L-4r + 4r-cl) (9) 

creates an electron in the d-channel. Notice that the 
spin exchange involves processes where the electrons "hop 
and flip" between neighboring orbitals. If large repulsive 
interactions are present in the conduction sea, then an 
electron can no longer "hop and flip" onto a site that is 
already occupied. This restriction means that creation 
operators must be replaced by Hubbard operators, 

c Jff — ■+ C ja (l - Tlj-a) = Xj>, (10) 

To see how this modifies the spin-exchange processes one 
can use a Gutzwiller approximation, 



where x is the concentration of carriers per site. This ap- 
proximation correctly describes the complete suppression 
of hop and flip processes in the limit where x = 1. With 
this replacement the transformed Kondo interaction de- 
velops three new scattering channels, 

Hi = JiS ■ (*V*«i) 

+ J 2 S • (*t s a* s + &p.** Pm + * f PB a* p J, (12) 

where J x = (1 - f ) J, J 2 = f J, 

*L = 2( c la+ c L+4<r + c L)> s-channel 

Vfrf - J_C r t _J \ ~\ 

lr : f / y ' p-" as) 

create electrons in the secondary channels. Electrons in 
the secondary channels are able to exchange spin with the 
local moment even though they do not hybridize with it. 

In more complex Uranium heavy fermion systems, 
intra-atomic interactions plcfth a vital role in opening 
up second-channel couplingsHd In Uranium atoms, the 
Hunds interactions have the effect of suppressing fluctu- 
ations in the "shape" of the localized orbital, so that 
electrons scattering off a localized orbital tend to ex- 
change spin, whilst preserving their orbital quantum 
numbers. In a tetragonal crystal for example where the 
low lyirufcstate of the f 2 ion is a magnetic non-Kramers 
doubletHil 

|±)=a|±l)+/3|T3>, (14) 

spin fluctuations within this doublet involve the exchange 
of spin with conduction electrons in two different "shape" 
channels, with equal Kondo coupling constants. 

III. TWO-CHANNEL KONDO LATTICE MODEL 

This discussion motivates us to examine how addi- 
tional spin exchange channels might modify the physics 
of a Kondo lattice. To this end, we shall consider a Kondo 
lattice model where two orthogonal scattering channels 
dominate the spin exchange process: 

H = ^ ^a^-kaCWa + ^ J T^T 3 Cr^ T j ■ Sj, (15) 

where . = (V4\fp ^Til)> = 1 ; 2) is a two component 
spinor 

tfrja = N7 h $rkC + k CT e- ikR >, (16) 

k 

that creates an electron at site j in one of two orthogonal 
Wannier states, with form-factor $rk- Here N s is the 
number of sites. We shall show that channel interference 
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Jj : s-channel 



J 2 : d x 2. y 2 -channel 




FIG. 3. Illustrating spin coupled to electrons via an s, and 



a d x 2_ y 2 channel. 



Jl : h yi -channel 



J : p - channel 

2 z 





FIG. 4. Illustrating spin coupled to electrons via a primary 
fxyz and a secondary p z channel. 



becomes strong when two channels have the same spatial 
parity. 

A simple example of our model is a two-dimensional 
tight-binding lattice of conduction electrons, where 



gk = —2t(cosk x + cosky) — /i, 



(17) 



and p, is the chemical potential, interacting with a local 
moment at each site in an s and a d-channel, so that 



$ik = 1, 



(s-channel) 



= (cos k x — cos k y ), (d x 2_ y 2 -channel). 



(18) 



as shown in Fig. ||[ 

A slightly more appropriate example would be a three- 
dimensional lattice, where 

ek = — 2t(cos k x + cos k y + cos k z ) — Mi (19) 

with a local moment at the center of each cube of atoms 
interacting in a primary f xyz channel, and a secondary 
p z channel: 



sin k x sin k y sin k z , (fxyz -channel) 
*2k = \/2smk z , 



y Z , 2 q. 

(p z -channel). 



as shown in Fig. |ij 



Unlike earlier treatments of two-channel Kondo prob- 
lems, our model involves a single conduction electron 
band, and there is no globally conserved "channel quan- 
tum number" . In a heavy fermion system, the orbital 
channels are locally well defined, but an electron scat- 
tering in one channel at one site, can then scatter in a 
different channel at a second site. This is important, for 
it can lead to interference effects between the Kondo ef- 
fect in different channels which are completely absent in 
models with an artificial channel quantum number con- 
servation. To illustrate this important point, we shall 
contrast the properties of our model with the channel 
symmetric "control model" 



H c = ekcVko-crkfr + Jr^rjcrcrj ■ Sj. 



klV 

where now 



c f r> 



= N s s cVfae" 



(21) 



(22) 



(r = 1, 2). In the control, electrons in different channels 
do not interfere, and we shall show that this prevents the 
development of composite pairing. 



IV. COMPOSITE PAIRING INSTABILITY OF 
THE ONE-CHANNEL KONDO LATTICE 

To examine the effect of second-channel couplings, we 
introduce the composite operator 



(23) 



This operator transfers singlets between channels by 
adding a triplet and flipping the local moment. To see 
how this works, consider a single site, where 



l*-a) 



1 

i 

V2 



^tIO-^iIT) 



(24) 



represent Kondo singlets in channels one, and two respec- 
tively. The action of the single site composite operator 
A^ = itp\aa2^ 2 • S is as follows: 



At|*. 1 ) = 2^ u V t iTl*-a>, 
At|* a2 )=2^t 2T ^t 2i | $al ) j 



(25) 



showing that the composite operator A^ transfers a 
Kondo singlet between channels, leaving an electron pair 
behind in the channel formerly occupied by a Kondo sin- 
glet. 

We now show how channel interference in the one band 
model causes the susceptibility of this composite opera- 
tor to develop a BCS-like divergence in the Fermi liquid 
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ground-state. Suppose that J 2 = and J\ is sufficiently 
large for a Kondo effect to develop in channel one. In the 
corresponding Fermi liquid ground-state |$), the com- 
posite pair susceptibility is given by 



XA 



E 



E<& — E\ 



)+(A#At)| 



(26) 



To evaluate the matrix elements appearing in this ex- 
pression, we need to decompose the composite opera- 
tor in terms of quasiparticle operators. The essence of 
the Kondo effect is the development of Fermionic bound- 
states between the local moments, and the conduction 
electrons. At low energies, the operator (Sytr a p)ipip then 
behaves as a single bound-state fermion, represented by 
the contraction 



(Sj ■ o a p)i>ip{j)= zf 3 , 



(27) 



where z is the amplitude for bound-state formation. By 
making this contraction, we imply that in all matrix el- 
ements between low- lying excitations \a) and \b) of the 
Fermi liquid, (Sj • (Taptyip can be replaced by a Fermi 
operator, as follows: 



(a\(Sj ■ a a p)*l>ip(j) \b) = z(a\f ja \b). 



(28) 



It is the contraction of the exchange term which gives 
rise to a resonant hybridization between f and conduction 
electrons 

J 1 [^ lj {S~7)yj lj +H.c] = Ji^ii/j + H.c.]. (29) 

so that at low energies, the Kondo Hamiltonian can be 
replaced by an effective Anderson model. 

The low energy eigenstates of the one channel Kondo 
lattice model are then an admixture of electron and com- 
posite fermion ako- = cos (5kCker + sin(5k/ker, with Hamil- 
tonian H* = 2kcr -E'ki^ko-akcr- The volume of the Fermi 
surface now-.oou.nts both the conduction and composite 
f-electrons.Em In the one band model, the conduction 
and composite f- electron share a single Fermi surface 
and they may be decomposed as follows 



C k(T = COS^kpflka 

/kcr = sin£ kp a k(T 



(30) 



where the high energy components that do not affect the 
low-energy matrix elements have been omitted . Near the 
Fermi surface the scattering is resonant and <5k F ~ 7r/2. 
Moreover, the small conduction electron admixture at the 
Fermi surface must reflect the symmetry of the screening 
channel, so that cos<5k F oc $ik F - 

We can now apply the contraction procedure to eval- 
uate the matrix elements of the composite operator. Let 
us begin with the control model. Applying the contrac- 
tion procedure we obtain 



(A|At|$) c = -iJ2(M Sj ■ (4j <ra 2 ct 2j )|$) 

j 



(31) 



k.cr 



respectively 



In the control model, c| k and /_ k ~ a)_ 
create light and heavy electrons on completely different 
Fermi surfaces. The mismatch between the volume and 
the dispersion of the Fermi surfaces for channel one and 
two assures that the excitation energy E\ — E$> = e^ + E^ 
is always finite: 



(A|At|$) c ex £ k( ^<A|ct 2k(T at_ k _ CT |<I>) 



E x -E^ = e k + £ k >0. 



(32) 



The channel susceptibility xa is consequently finite. We 
conclude that with perfect channel symmetry, a small 
second-channel coupling is irrelevant. 

Now let us remove the channel symmetry and return 
to the physical model. Now we have 



<A|At|$)=z^(A| S r (^aa 2 ^ 2j )\<£) 

j 

= zJ2^(X\crcij^ k . a \^. (33) 



Unlike the previous case, this pair creation operator can 
be decomposed in terms of quasiparticles on a single 
heavy Fermi surface. Transforming to quasiparticle oper- 
ators using eq. (|30| ) introduces a factor cos((5k) sin(<5k) ~ 
$ik into the sum, so that 



(A|At|$) oc E kCT ff $ i-k*2k(A|at kCTa t_ k _ CT |$) 



E\ — E<s> — 2Ek 



(34) 



This relation describes the decomposition of the com- 
posite pair operator in terms of the low-lying quasipar- 
ticles. Notice that the matrix element is proportional 
to $i_k$2k, showing that this amplitude involves an in- 
terference between the two channels. Furthermore, the 
two form-factors must must have the same parity, or the 
composite operator vanishes on the Fermi surface. Since 
the excitation energy, 2E^ vanishes on the heavy Fermi 
surface, it follows that there are now a large number of 
zero modes for the transfer of singlets between channels. 

It follows that the composite pair susceptibility xa now 
contains a singular term. Substituting the above results 
into the general expression for the composite pair suscep- 
tibility, we find 



XA oc 



($ik$ 



2k) 



k 



2E U 



00, 



(35) 
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AI«D> = 



ch. 2 



A|<D> = 




FREE MOMENT 



AE = 2E k 



AE=8 k+Ek 




ch. 1. 

FIG. 5. Action of composite operator on heavy Fermi liquid 
creates: (I) a pair of heavy fermions (channel interference ) 
and (II) a heavy and light electron (channel conservation). 



which diverges logarithmically in the thermodynamic 
limit. We see that once channel symmetry is broken, the 
composite pair susceptibility xa is directly proportional 
to the BCS pair susceptibility of the heavy quasiparti- 
cles, where the symmetry of the channel is given by the 
product of the two screening channels. 

This has immediate consequences for the effect of a fi- 
nite J2 on the Fermi liquid ground state. Once channel 
symmetry is broken, the susceptibility to transfer singlets 
by creating composite pairs diverges. Any finite J 2 will 
polarize the transfer of singlets into channel two, thereby 
coupling J2 to this divergent susceptibility. Thus the loss 
of channel symmetry causes a coupling to a second chan- 
nel to become a relevant perturbation. This will force 
J2 to scale to strong-coupling. A similar conclusion will 
hold when J2 is large, and J\ is small. The simplest 
way to connect up the renormalization flows in the vicin- 
ity of the strong-coupling Fermi liquid fixed points with 
the flow away from the weak coupling fixed point is by 
hypothesizing the presence of a new attractive Kondo lat- 
tice fixed point that is common to both channels. (Fig. 



V. SU(2) FORMALISM 

The key to the development of a field theory for com- 
posite pairing, lies in the use of the Abrikosov pseudo- 
fermion representation for the local moments 



Sj — f ja 
n f (j) = 1. 



v 2 , 



a/3 



fj0> 



(36) 



J r J 2 =oo 



FIG. 6. Conjectured renormalization group flows for 
the co-operative two-channel Kondo effect. The Fermi liq- 
uid formed in channel one or two is unstable to common a 
two-channel state with composite order. 



Since f-charge fluctuations have been removed, the KLM 
is defined within the sub-space constrained by the 
(Gutzwiller) requirement rif = 1 at each site. The ab- 
sence of f-charge fluctuations is manifested as a local 
Slj^j!) gauge invariance of the Heisenberg spin operator 



J jc 



pi<P ft . 

cosQp j a + sgnasinOfj-o-. 



(37) 



To illustrate this feature, consider the spin raising oper- 
ation S+. This process can proceed by first annihilating 
a down electron, then creating an up electron, written 
Sj~ = p jifji- Alternatively, it can proceed by first cre- 
ating an up electron, forming the n/ = 2 state, then 
annihilating a down electron, written Sj~ = — • 
In fact, one can accomplish the spin raising operation by 
an arbitrary linear combination of the above: 

S+ = (cos 6p n + sinOftt ) (cos Qfa - sin6p n ) . (38) 

In other words, there is no distinction between a particle 
or a hole when all charge fluctuations are removed .0 

The SU(2) symmetry implies that the constraint n/ = 
1 is actually component of a triplet of local "Gutzwiller 
constraints" 



= 0, 



which can be written in the compact form 

/ V£ = 0, 



(39) 



(40) 



G 



where f^j — (fj\,fji) is a Nambu spinor, and r = 
T2, T3) represents the triplet of Pauli matrices. The 
hrst two constraints are particularly important in any 
consideration of a paired state, providing the main driv- 
ing force for anisotropic pairing. 

The partition function for our model is given by Z = 
TrlPce^ 13 "] where P G = Wjin^-n^) 2 is the Gutzwiller 
projection for one f-spin per site. Following earlier work, 
we rewrite the Gutzwiller projection as an integral over 
the SU(2) group 



d[Wj] g h 



(41) 



where §j = i w ih is the SU(2) operator, Wj — 9jhj-T, 
( 9j e [QJtt]) and d[W] = sin 2 6d6dh/ {An 2 ) is the Haar 
measures over the SU(2) group. Introducing this into 
the partition function permits us to write it as a path 
integral 



Z = 



V 



r, rxri - {"{Ci+H)dT 

[j, c, W\e Jo 



(42) 



where 



£l = E J&Ck + E PiVr - W i)fi> ( 43 ) 



is the Berry phase. 

The antiferromagnetic interaction between the local- 
ized moments and the conduction electrons can be de- 
coupled in the particle-hole channel, as follows 



Jr [<rr ■ S - i] = - y W, a r }, 



where 



a r 



E/ f ^ra 



(44) 



(45) 



The SU (2) gauge symmetry guarantees that there is in 
fact, a continuous family of ways to decouple the interac- 
tion. Thus, by making the transformation f a — > af- a \ 
we can decouple the interaction in the Cooper channel, 
as follows, 



where 



J r [<r r .S-i] =-±{brKb r }, 



br = E a f-^o 



(46) 



(47) 



We now decouple the interaction simultaneously in both 
channels, by first writing 



Hi = -'^-[{a r \a r } + {b r \b r }], 
then decoupling each term as follows 



(48) 



- ^W, ar} - [ar f ^ r + h.c] + ^-V r *V r , 
4 Jp 

-4-{°r f ,M^ [6r t A r +h.c] + -|-A r *A r . (49) 
4 Jp 

It is convenient at this point, to introduce a Nambu 
spinor representation for the conduction electrons 



Ck = t CkT 



(50) 



The corresponding spinor for the localized electron Wan- 
nier states is 



ikR, 



(51) 



where p is the parity of the form-factor <I>rk = P'I'r-k- 
The decoupled interaction can now be written in the sym- 
metric form (Appendix A ). 



1 



J r (S ■ a r ) - [fVrVT + H.c] + — Ti[V\Vr], (52) 

where Vr is directly proportional to an SU(2) matrix gr 

r 



V r = 



V A 

A* —V* 



= iV L 9v, 



The integration measure for Vr is 

d[V T ] = dV r dVrdA r dA* r . 



(53) 



(54) 



Repeating this decoupling procedure at each site in the 
path integral, enables us to write 

Z = J £>[/, c; W, V]e~So ^+ H ) d ^ 

H = E e £ c V 3C fc + /// ' 



Hj = Ejl/^Vr^r, + H.c] + j-Tr[V^ n V n ]j. 



(55) 



VI. GAUGE FIXING 

Our model now has the following time dependent 
SU(2) gauge invariance 



fj -> 9jfj, 
V T] -> gjV T j, 



W 3 ^g {W 3 +id T )g^ 



(56) 



associated with the absence of f-charge fluctuations. 

When we develop a saddle-point expansion for the 
functional integral, we need to deal with the local zero 
modes associated with this gauge invariance. Following 
standard gauge theory practice, this means that we need 
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to fix the gauge. We choose the "radial gauge" , where the 
Kondo matrix V is proportional to a unit matrix in the 
channel with the largest Kondo coupling constant. This 
is the SU(2) analog of the radial gauge used by ReacLandi 
Newns in their U(l) treatment of the single channelEIl'Eij 
Kondo lattice. Suppose that 



Vi(i) = ivi(j)h 



(57) 



where Vi(J) is real and hj is an SU(2) matrix. To fix the 
gauge, we absorb hj into a redefinition of the fields by 
setting gj = fv j(r) and making the gauge transformation 
(pq). In the radial gauge, 



Vij(t) = wy(r)l, 

Wj(r) = h^{r) {Wj + id T ) hj(T), 



(58) 



so the formerly static field Wj is elevated to the status 
of a dynamic field Wj(r). The measure for the bosonic 
fields inside the path integral is now 



d[W,V] = M 3 dM[V 2 ]d 3 W, 
at each site and time-slice. 



(59) 



VII. LINK BETWEEN COMPOSITE ORDER 
AND CHANNEL INTERFERENCE 

Under the local SU(2) gauge transformation fj — > gjfj, 
( [56| ) the two amplitudes Vy and V2j which describe the 
Kondo effect at site j transform in precisely the same 
way. The only gauge invariant term we can form from a 
single channel is trivially proportional to the unit matrix: 
VljV V j oc= 1. But in a two-channel Kondo problem, the 
interference term 



2j rij 



(60) 



is also gauge invariant, since V^2jVij — ► V'ij9^j9jV%j = 
V'2jVij and for this reason, is expected to have a sim- 
ple physical significance. To identify the meaning of the 
interference term we introduce a source term into the 
Hamiltonian that couples to it 



H^H + ^Tr[V 2 /V : 



h.c] 



(61) 



and 



Fj = it\jOi) 2 j ■ Sj, 
Aj = ijjijia y a^2j ■ Sj 



(64) 



represents the composite order in the particle-hole, and 
particle-particle channels respectively. By comparing 
© and (||), we obtain a special relationship between 
the inter-channel interference and the composite order, 



V^jVrj = - 



J1J2 



F A 

-At Ft 



(65) 



Notice incidentally that the off-diagonal terms are odd 
under interchange of the channel index. 

We thus learn that if the Kondo effect develops coher- 
ently in two channels, composite order develops. This 
enables us to understand why composite order develops 
critical-correlations in the symmetric two-channel Kondo 
model.U In a lattice, true long-range order becomes pos- 
sible. 

Let us briefly consider the possible phases that might 
develop. If Vr develops a finite amplitude in both chan- 
nels then the composite order in the ground-state will 
have the form 



MF|V> (V>|A|V> 



J1J2 



V\,V :r (66) 



Suppose the amplitudes of Vr are constant, then in the 
"radial gauge" 



V 



ivil, 



V 2] = iv2e~ 



(67) 



where the vector rij develops a vacuum expectation value. 
The composite order matrix is then 



wm) mam 

-<V|At|^) (V>|Ft|^> 



= .\/, : ('" " T . 



(68) 



where M Q = 2v\V2/{J\ ^2)- Two kinds of phase are pos- 
sible: 

• Composite magnetism where rij = z. In this 
phase, the order parameter matrix is diagonal, and 

(iP\^ lj <rf 2 j-Sj\iP)=M e i ^. (69) 



where aj = a® + iotj ■ t is a unitary matrix, with four 
real coefficients (a , a) at each site . If we now reverse 
the Hubbard Stratonovich transformation, by integrating 
over the fields Vr (Appendix B) the Hamiltonian acquires 
the additional term 



H ^ H + Tr[Mj^aj + h.c], 



where now 



A^t 



J1J2 



F A 
-At Ft 



(62) 



(63) 



This phase breaks time-reversal symmetry, forming 
an orbital magnet where the spin becomes corre- 
lated with electrons in two orbitals. 

• Composite Singlet pairing where 4>j = ir/2. If 
n(x) = cos9(x)tc + sin0(x)y, whereupon 



A s (xj) = i(il>\Tjjij(jyCnp2j ■ Sj\tp) = iM a e 



- io , 



(70) 



The second possibility is particularly interesting, because 
the composite pair susceptibility diverges in the Fermi 
liquid phase. This is the main topic of the of the paper. 
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VIII. MEAN FIELD THEORY OF THE 
COMPOSITE PAIRED STATE 

We now develop a mean-field theory for the uniform 
composite paired state. With this theory, we show that 
the two strong-coupling Fermi liquid phases of our two- 
channel Kondo model share a common instability into a 
phase with uniform composite order. 

We seek a uniform solution, where all mean-field pa- 
rameters have no dependence on position. In this case the 
mean-field Hamiltonian is most compactly represented in 
momentum space, as 



H 



AIF 



£kT 3 V f k 
V k W T 



where now 



V k = V^ik + V 2 * 



2k- 



(71) 



(72) 



Strictly speaking, here we should have written the form 
factors as $rkT 3 , to take account of the possibility of an 
odd-parity scattering channel. However, provided both 
channels have the same parity, we can always cast Vk in 
the above form. For even parity channels, $rkr 3 = ^rk 
directly. For odd-parity channels, 'I'rkTa = ^rkTs, but in 
this case the t 3 can be absorbed by a gauge transforma- 
tion fj -> T 3 fj, 

To examine uniform pairing, we shall take 



Vi = it'll, 
V 2 = v 2 n ■ t, 
W = Ar 3 , 



(73) 



where n = cos6>y — sin#x describes the phase of the 
composite pairing. For convenience we shall take n = y, 
so that 



V k = ivik + v 2k r 2 , 



(74) 



where we have introduced the notation t>rk = ^r^rk- 

Ostensibly, our mean-field theory is that of a BCS su- 
perconductor, with Hamiltonian described by 



W(k) = 



£kT 3 

v k 



vt k 

W T 



(75) 



However, there is one important distinction: here the 
pairing takes place between charged conduction electrons 
and the neutral f-spins, and is merely a manifestation of 
the formation of composite pairs. For this reason, it is 
actually not possible to say whether the pairing is channel 
one, or in channel two. In the gauge we have chosen, the 
scattering in channel one is "normal" and pairing takes 
place in channel two. But suppose we make the gauge 
transformation ( |56|) with gj = — ir 2 , then 

V k = iviu + V 2k T 2 iT 2 Vk = Vl^T 2 - iv 2 k, 



W = Ar 3 — > iT 2 W(— it 2 ) 



-Ar 3 , 



(76) 



which transforms the Hamiltonian to one which is now 
pairing in channel one, and "normal" in channel two. We 
are forced to recognize the superconductivity can not be 
identified with either channel, but instead derives from a 
coherence between the two channels. 

Suppose we now integrate out the f-electrons: now we 
find that the conduction electron Green function has the 
form 

GOsJ-^w-ewTa -£(«), (77) 

where k = (k, ui) and where the self-energy term 

E(/s) = V\{u; - Ar 3 )- 1 V k (78) 

describes the resonant scattering off the quenched local 
moments. If we expand the self-energy, we see that it 
contains both normal and anomalous components 



(79) 



Notice that although this self-energy contains off- 
diagonal terms, it is invariant under the local SU(2) 
gauge transformations. The normal components are 
channel-diagonal 



Sat(k) 



'ik 



y 2k 



(80) 



u) — At 3 u) + At 3 ' 
but the anomalous terms depend on channel interference 



2Avik«2k 



(81) 



and are directly proportional to the composite order pa- 
rameter A. For A /: 0, composite order induces con- 
ventional pairing amongst the conduction electrons. We 
shall later see that the independent existence of the com- 
posite order means that a finite Meissner stiffness devel- 
ops even when A = 0, and conventional pairing is absent. 

The mean-field free energy per site can be written as 
follows 

F = -T/N. J2 Tr H^n - W(k)] + 2 J- (82) 
kiw^ r 

The eigenvalues of the mean- field Hamiltonian H(k) oc- 
cur in two pairs (— utk v , w^), where rj = ±, corresponding 
to two bands of quasiparticle excitations. We may rewrite 
the characteristic determinant of Tt(k), Detfo; — W(k)], 
in terms of the G(k) 



Det[w - H(k)} = Det[G(«) _1 ](w a - A 2 ). 



(83) 



To evaluate the determinant, we write G(k) 1 in the fol- 
lowing form 



G(k)- 1 = (A-Btz + Ct^-X 2 )- 1 , 



(84) 
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where 



A = uj{uo 2 - A 2 - vl), 
B = e k (u; 2 -X 2 ) + \v 2 _, 
C = 2Awi k W2k, 



(85) 



uf k , and v k _ = v 2 k ± vf k . The eigenvalue 



equation Det[w — W(k)] = then becomes 
(A 2 — B 2 — C 2 )(u 2 — A 2 ) -1 = 0. 
Expanding this expression, we obtain 

Det[w - W(k)] = - 2w 2 a k + 7 £, 

where 



(86) 
(87) 



and 



«k = «k + 5(4 + A 2 ), 



7k = \l (Ae k - t>£_) 2 + (2vi k i;2k) s 



The eigenvalues of 7i(k) are thus given by 



w k± = ^/a k ±(a k - 7 D 1/2 
and the mean-field free energy is 
2T 



k.r/ 



2 cosh( / 5^k?7/2) 



r=i,2 



Jr 



(89) 



(90) 



(91) 



By minimizing the Free energy with respect to A, v% and 
V2, we can now determine the mean-field phase diagram. 



IX. PHASE DIAGRAM 

We now discuss the mean-field phase diagram. There 
are three types of stable mean-field solution: 

• Normal phase. Ui or V2 7^ 0, V1V2 = 0. At 
high temperatures, either V2 or v% is finite finite, 
signalling a Kondo effect in the stronger channel. 
There are thus two types of normal phase with dif- 
ferent Fermi surface geometry, depending on which 
channel is the strongest. Suppose v 2 — 0, then 
the normal state spectrum ( P0[) attains the simpler 
form 

% — ► £k± = - [(e k + A) ± \J (e k - A) 2 + 4v 2 k ] , 



corresponding to a band formed by an admixture 
between the conduction electrons, and the compos- 
ite f-electrons in channel one. This phase describes 
a heavy fermion metal. 



• Gapless composite paired state. v%V2 7^ 0. 
In the generic composite paired ground-state, the 
quasiparticle excitation contains nodes. The con- 
dition for gapless excitations is 



7k = 0, 



which implies that 

(Ae k - v 

Ulk«2k 



k-) = 0, 
= 0. 



(92) 



(93) 



The first condition defines the locus of points on 
the underlying Fermi surface, whilst the second- 
condition defines the nodes of the order parameter. 
Gapless quasiparticles form at the intersection of 
the order parameter nodes with the Fermi surface. 
This occurs when one, or the other channel is dom- 
inant and in this case, the conduction propagator 
(\Fi\) can be written 



G{k) 



Z^w-Efo- A k n], 



(94) 



where the gap symmetry is determined by the prod- 
uct of form factors, 



A k = A $i k $2k, 



and 



(95) 



(96) 



is a mass-renormalization constant, and E k = 
Z k (e k — w k _/A), A D = 2Z k viv 2 \~ 1 describe the 
kinetic and pairing contributions to the quasipar- 
ticle energy. The heavy electrons in this state are 
paired, and thus correspond to an anisotropic BCS 
superconductor a spectrum w k = (E k ) 2 + (A k ) 2 . 

• Gapped composite paired state. V1V2 7^ 0. 
The gapped composite paired phase occurs when 
( = J2/J1 lies between two critical values, 



Ci < J2/J1 < C2- 



(97) 



If the weaker coupling constants is increased to 
the point where it is comparable with the stronger 
channel, the underlying Fermi surface collapses 
around the zone center, causing the the nodes to 
mutually annihilate. At a a still larger coupling 
constant, the nodes reappear at the zone corners. ( 
Fig.^) This phenomenon is perhaps easiest to visu- 
alize when the conduction sea is half-filled. In this 
case, the normal-state is a Kondo insulator with 
no Fermi surface.EJ The mean-field theory predicts 
that when J2 exceeds a critical value, composite 
pairing can take place forming a pure composite 
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FIG. 7. Evolution of the gap nodes(open circles) and the underlying Fermi surface in the composite superconductor, as the 
ratio £ = J2/J1 evolves from zero, to infinity. In the intermediate region, where £1 < J2/J1 < (,2, the underlying Fermi surface 
collapses around the zone center, and the nodes annihilate one-another to produce a gapped state. 



paired state. Although the quasiparticle spectrum 
is gapped, as in a Kondo insulator 



4< k 



±ck], 



(98) 



the composite order parameter A = 2v\V2/ (J1J2) 
is finite and there is a superconducting response. 
When this gapped state is doped, it preserves its 
gap. 

In the paired phase, the mean-field equations are given 
by the three conditions 



dF 
0A' 



OF 
dvi ' 



OF 
dv 2 



= 0. 



(99) 



The first of these equations imposes the constraint 
/T3/ = 0, whereas the second and third determine the 
magnitude of the Kondo effect in the two channels. Writ- 
ten out explicitly, the mean-field equations are 












1- 


3) 











A = 



'akA - e k (e k A - v 2 . 
l(e k + A) 2 $? k 



(100) 



i(e k -A) 2 $ 



2k 



Suppose channel one is dominant, then in the normal 
state vi = 0, which yields two equations for the normal 
state 



E th 

k,n=± 



e k - A 



2 



0. 



= X k(T), (101) 



where 



1 



k,r)=± 



th 



• {3E kr] N 



*?k 



E, 



Eu- 



k— rj 



(102) 



By setting v 2 = + in the full set of mean-field equa- 
tions, we find that the transition temperature into the 
composite paired state is given by 
2 

Xc(T c ), 



where 



and 



Xc(T) 



1 

~Ns 



E th( 

k,7?=± 



A 2 _ ({,2 

yi k?j — v 2k 



1 + 



(e kF -A) 2 



- C/ k4 



(103) 



(104) 



(105) 



ement A klJ 



identified as the matrix el- 
\(k F t,—kp ||A|$)| 2 associated with the 
action of the composite pair operator on the Fermi liquid 
ground-state |$) . In all the above expressions, we have 
simplified the algebra using the identities 



Ekn - E k - V = t]\J (ek - A) 2 
Ek n + Ek- V = e k + A. 



(2wik) 



(106) 



We now discuss the detailed phase diagram that results 
from the mean-field equations. 
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A. Instability of the heavy electron metal 

We have argued that the presence of a heavy fcrmi 
surface leads to new zero modes for the transfer of sin- 
glets between different screening channels, so that when 
a second channel develops a finite coupling, a composite 
pair instability immediately results. We now describe in 
detail, how this result emerges naturally from our mean- 
field theory. 

Let us begin by setting the scale of the Kondo temper- 
ature within this mean-field theory In mean-field theory, 
the cross-over into the Fermi-liquid regime is crudely de- 
lineated by a mean-field phase transition. We use this 
temperature as a definition of the single site Kondo tem- 
perature. At the "transition temperature", v\ = + , so 



Xk(T 



K 



N, 



1 



(107) 



where we have replaced the momentum sum by an in- 
tegral over energy, so that N(0) is the density of states 
at the Fermi surface, and ($f k ) denotes a Fermi sur- 
face average of the Form factor. With this definition the 
mean-field Kondo temperature takes the form 



'-K1 



£) e -l/JlJV(0)<*? k >_ 



(108) 



This quantity sets the characteristic size of the mean-field 
parameters in the normal state: 



A ~ (Vi) 2 N(0) <~ Tki- 



(109) 



The composite pair instability of the normal state is 
directly related to a divergence in the fluctuations asso- 
ciated with the Kondo effect in channel two. To see this, 
we expand the mean-field expression for the Free energy 
to Gaussian order in V 2 = v 2x t x + v 2yT y , which gives 



F = F + \v 2 \ 



T 2 - Xc 



(110) 



where v 2 = v 2x + iv 2y and \c is given above. From 
this result, we can read off the fluctuations in the order 
parameter, 



(5v 2 6v2) 



(111) 



But since the composite order parameter is given by 
SA = 2v\5v 2 /(JiJ 2 ), it follows that the composite pair 
susceptibility is given by 



X a = [(SA6A*) - (5A5A*) XC=0 ]/T 
vi\ Xc 
JiJ l-(J 2 Xc/2)' 



(112) 



The denominator in this expression vanishes at T c , ex- 
plicitly confirming that the composite pair susceptibility 



diverges at the mean-field transition between the normal, 
and the paired state. 

To gain some insight into the composite pair instabil- 
ity, we divide the bare composite susceptibility, xc into 
a "high" and a "low" energy component 

Xc= J2 + J2 {■■■} = Xh + Xl, (113) 

where the former describes the local Kondo effect in 
the weaker channel, the latter, the channel interference 
taking place on the heavy Fcrmi surface. At energies 



EM 



2K 



(114) 



Suppose the heavy Fermi surface lies in the lower band, 
on the Fermi surface, 



J ik 



so that 



4A 2 ^($ lkF $ 2kF )2 



A 2 

A k FJ; - 



(A 2 +< F ) 2 



(115) 



(116) 



When we replace the momentum sums by energy inte- 
grals, we must remember that the density of quasiparticle 
states is enhanced by a factor 



A^ F (0) = iV(0) 



<ie k 



dE u 



(A 2 



k=k« 



A 2 



^ k) A(0). (117) 



The energy scale separating the two regimes is the Kondo 
temperature for channel one, Tki- With these results, 
approximate expressions for the high and low energy con- 
tributions to the composite susceptibility are 



Xh 



2N(0) 
Xi 



JT K1 e 



( TK1 ^ / ( V l*lk*2k) 2 

2N{0) J T E 

[ Tki dE _ 2 . 



(A 2 +« 2 k ) 



(118) 



The expression for xi was simplified by noting that the 
dominant contribution to the second term occurs in the 
regions far from the node in the order parameter, where 
vik » A. The sum of the two expressions then yields 



Xc « 2N(0) 



ln(_^) + 21n(^i) 



(*lk>- (H9) 



The Gaussian coefficient of v 2 in the Free energy is thus 
given by 
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-r ~ Xc 



2iV(0)($ 



2k) 



92 K T K1 - 



21n 



■ Tki 
■ T 



(120) 



where 92 = iV(0) J2($2k) IS t ne dimensionless Kondo cou- 
pling constant for channel two. The first logarithm in this 
expression describes the renormalization of the coupling 
constant in channel two down to the energy scale Tki'- 



92(T 



Kl 



92 y T K1 h 



(121) 



the second logarithm describes the subsequent renormal- 
ization of 91 at temperatures below Tki- In a two- 
channel single impurity model, the second logarithm 
would be entirely entirely absent because the Kondo ef- 
fect in channel one cuts off any further renormalization 
in channel two. Here we see that the constructive inter- 
ference between the two channels in the lattice actually 
over-compensates for the Kondo effect in channel one, 
producing a logarithmic renormalization at low temper- 
atures which is twice as large as at high temperatures. 
The co-operative Kondo effect thus develops at a tem- 
perature which is higher than the Kondo temperature for 
an isolated channel two. We may rewrite the Gaussian 
coefficient as 



2 

J2 



Xc 



2 

T 2 



2N(0) 



<*2k> 



-2JV(0) 



In 



Tk2- 



M 



4JV(0)($i k )ln(- 



Tki 
T 

). 



\JTk\Tk2 

where we have used the definition 

T K 2 = De- 1/{N(0)M ^ ) \ 



(122) 



(123) 



to absorb the coupling constant Ji- In this rough ap- 
proximation, the composite pairing instability occurs at 
a temperature 



T c ~ \JTk\Tk2- 



(124) 



When J2/J1 « 1, this same scale sets the size of A Q . 
Fig. |§| illustrates the phase diagram calculated numeri- 
cally for a two channel Kondo lattice with "s" and "d- 
wave" screening channels: 



$ik 

$2k 



1, 

(cos k x — cos k y 



(s-channel) 

{d x 2 _ y 2 -channel) . 



(125) 



as shown in Fig. |^. As expected, the composite pair 
instability occurs at the highest temperature when the 
two channels are most evenly matched. 



— — — — — — — jr — 


— — — 


■ FL 1 / \ 


FL 2- 


/ At^O n 




0.5 1 


1.5 


J a /J, 





FIG. 8. Phase diagram for a two channel Kondo lattice 
with "s" and "d-wave" screening channels. In the shaded re- 
gion, a co-operative Kondo effect in both channels gives rise to 
composite pairing and a quasiparticle gap with sxd — d-wave 
symmetry. 



B. Composite pair instability of the Kondo insulator 

Since composite pairing is an intrinsically local process, 
the presence of a heavy fermi surface is not a necessary 
requirement for the formation of the paired state, but in 
its absence, the instability requires the second-channel 
coupling constant to exceeds a critical value. 

This is precisely what happens where the conduction 
band is half filled, for in this case, the normal state of 
the Kondo lattice is a "Kondo insulator" , with no Fermi 
surface and a gap to charge excitations. The Kondo in- 
sulator is particle-hole symmetric, so in the mean-field 
theory, A = 0, so that the excitation spectrum simplifies 
to the following form 



/EkN . 2 , £k 



(126) 



By minimizing the mean-field Free-energy with respect 
to variations in v\ and v%, and setting T = 0, we obtain 
two mean-field equations for the ground-state 



d d k 



Hi 



Jx J (2ir) d ^2 + ( 2 « k ) 2 ' 

1 r d d k $ k2 

T 2 'J (27T)* ^ £k 2 i (2^' 



(127) 



(where d is the dimensionality). This composite paired 
phase will only be stable within a range of J2, J| < 
J2 < J|*. By setting V2 = + , we obtain two parametric 
equations for the second-order phase boundary between 
the Kondo insulator and the composite paired state: 
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1 
1 

J? 



d d k 



(2«i$i k ) 2 



H2 



( 27r ) d Ve k 2 + (2 Wl $ lk ) 2 ' 



(128) 



Beyond the critical value Jj**, the Kondo effect is no 
longer operative in channel two. By setting v% = + , 
we obtain two parametric equations for the second phase 
boundary. 



1 

l 

J 2 



d d k 



d d k 



(277 2 $ 2k ) 2 



* k2 



(27r) d v/ ek 2 + (2t, 2 $ 2k ) 2 



(129) 



We have calculated the phase diagram for a two- 
dimensional Kondo insulator in two dimensions, with 
dispersion e k = — 2t(c x + c„) and a Kondo cou pling in 



the s- and d-channel, as shown in equation (125). When 
•h < J% , the Kondo effect in the s-channel leads to Kondo 
insulator. By contrast, when J\ << J 2 , a "nodal semi- 
metal" forms in channel two, with a gap which vanishes 
along the d-wave nodes of the form- factor $ 2k . Fi g. O 
show s the phase diagram obtained from equations (|12q ) 
and ( |129| ). Notice that the range of J 2 over which the 
composite paired state is stable, is negligible for small J\ 
and J 2 , but grows substantially as both J\ becomes large 
compared with t. 

This composite paired state is interesting, for its quasi- 
particle spectrum is essentially identical to the Kondo in- 
sulator, and furthermore, since A = 0, there is no anoma- 
lous component to the conduction electron self energy: 
there are no paired conduction electrons in the ground- 
state. Yet despite these similarities, the presence of com- 
posite order 



A(x) 



J1J2 



(130) 



means that this state is a superconductor, with a finite 
charge susceptibility. The calculation of the charge sus- 
ceptibility needs to be carried out subject to the con- 
straint. Expanding the Free energy energy to quadratic 
order in changes in A and the chemical potential, /x we 
have 



F n 



1 



[ Xw (<5/i) 2 + 2 X ^5fiSX + xxx(SX) 2 ] ■ (131) 



The constraint dF/ dX = implies that 



6X = 

so that we may write 
F = F D 

where 



X\x 



Xc , 



-Sfj,, 



(132) 



(133) 
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FIG. 9. (a) Phase diagram for two-channel Kondo insula- 
tor. "K.I " denotes the "Kondo insulating" phases. In the 
intermediate gapless phase both channels participate coher- 
ently in the composite pairing process, (b) Showing phase 
stiffness in the composite paired phase for the case Ji/4t — 1. 



Xc = Xfifj, 



Xlx 
X\\ ' 



(134) 



A rather laborious calculation (Appendix D) gives 



4w 2]3/ 2 » 



E 



k [(e k ) 2 +K] 3 / 2 ' 



XAA — XfifJ, + Xb, 



(135) 



where 



Xb 



^ WlkW2k ) 2 (5 + 4)- (136) 



Since \x»\\ < Xw and xx\ > Xc is positive, pro- 
vided both V\ and i> 2 are finite. In the next section we 
also confirm that the composite paired state has a finite 
superfluid phase stiffness, given by 



Ps = 



-y 



{viv 2 ) 2 ($l k V$ 2k - $ 2k V$i k ) 2 



(2« k+ ) 2 



(137) 



in the ground-state. 
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C. Gapped Composite Paired State 

In our simple two dimensional example, t>i k is never 
zero, so the condition ( |93| ) for gapless behavior become 

«?.. ^ 

(138) 



The modulus |(«i k /e k )| is smallest at the band-edges, 
where where |e k | = W± = 4i =F So if 




"lk 
W- 



< A < 



W+ 



(139) 



the paired state becomes gapless. Although a small resid- 
ual Cooper pair density is present, the state is essentially 
a pure condensate of composite pairs. To understand 
what happens at the critical values of £, it is instructive 
to examine the underlying Fermi surface, defined by the 
locus of points 



'2k 



A 



(140) 



As one approaches the critical value £ = £i, this Fermi 
surface collapses around the zone-center, forcing the 
nodes to mutally annihilate. In going from the region 
J2/J1 << 1 to the region J2/J1 » 1, A changes sign, 
and so that the system must always pass through this 
region of gapless composite pairing. In the center of this 
region, where A = 0, the state is a pure composite paired 
state. Once ( exceeds the value the nodes are re- 
born at the zone corners and the Fermi-surface re-appears 
along the zone-boundary (Fig. 0) 

To illustrate these conclusions, we have used the mean- 
field equations to calculate numerically the two lines 
where the gap vanishes. In Fig. [l], we summarize the 
results of these calculations, in a diagram where we have 
kept max(Ji, J 2 ) = 4t, and varied the ratio Jij J\- 



X. SUPERFLUID DENSITY OF THE 
COMPOSITE PAIRED STATE 

To confirm that the composite paired state is super- 
conducting, we need to compute the superfluid density 
p s . In the London gauge V • A = 0, the supercurrent is 
given by 



where 



j s = -QA, 



Qab — e 2 [Ps]ab — 



8 2 F 



<9A a <9A b 



(141) 



(142) 



In the presence of an electro-magnetic field, the electron 
kinetic energy and form factors acquire a dependence on 



the vector potential. Using a Nambu notation, we may 
write 



£k 
$rk 



^k— eAr3 ) 
^rk-eArs, 



so that the hybridization acquires the form 

^2T 2 $2k-eAr 3 -- 



V£ = Wl$lk-eAr 3 



(143) 



(144) 



The appearance of the vector potential in the form factor 
reflects the fact that the hop and flip motion of electrons 
around a local moment leads to current flow. In the 
London-gauge (V • A = 0) we may calculate by compute 
the second-derivative of the Free energy at fixed values 
of v\, V2 and A,cJ so the the only important part of the 
Free energy is the electronic component: 



F e = -T^Trln[zca„ - W A (k)]. 



(145) 



To second-order in the vector potential, we may write 

H A (k) = W(k) - J a A a + ^A a A b W 2 ab H(k) + 0(A 3 ), (146) 

so expanding the Free energy to second-order in A, we 
obtain 



F = F n 



K 

where 

g K = [iu n - ft(k)]- 1 = 

is the matrix propagator and 

J K = -v A n A {k) = 



G 

G fc 



V k £k 

V k V k r 3 



,G cf 
G ff 



raVkVn 




(147) 



(148) 



(149) 



i s th e current operator. The first and second terms in 
(147) correspond to the paramagnetic and diamagnetic 
components of the stiffness, as in a conventional super- 
conductor. Notice however that the current contains 
anomalous off-diagonal contributions that do not com- 
mute with the charge operator T3. Unlike a pure BCS 
superconductor, here the presence of composite pairs af- 
fects the current operator. 

The diamagnetic contribution to Q ab can integrated by 
parts, to obtain 



£ Tr[g K Vl b H K ] = - Tr[S K j^], (150) 

K K 

where 

J K = eV k H(k) = e Q , (151) 
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so that the full expression for the superfluid stiffness ten- 
sor is 

Qab = tYsI^GkJZGkJ*] - Tr[g K f K g K i b K }\. (152) 

K 

To gain some insight into this equation, it is instructive 
to evaluate the stiffness for the case of pure composite 
pairing at half-filling. In this case, the conduction elec- 
tron propagator commutes with the charge operator T3, 
so contributions to the stiffness which involve the con- 
duction electron component of Q identically vanish. For 
example, the cross-term between Vek and VV is 

Tr [G cc V a eG c/ V 6 Vr 3 - G cc V a er 3 G c/ V 6 V] = 0. (153) 

The only surviving terms represent the composite pair 
stiffness, which can be written in the form 



Qab - ^Tr [[G e/ VaV, r 3 ] [G cJ V b V, r 3 ] + H.c] . (154) 



From the Dyson equation, 
G(k) = G (k) + Q{k) 



V* 
V 



Go{k), 



(155) 



where Q is the propagator in the absence of any hy- 
bridization, it follows that 



G cf ( K ) = G(^)V t 
so that at half-filling (A = 0), 

G cf (K 

and hence 
[G c/ VV,r 3 ] 



1 



LO — At 3 



lo(uj - e k r 3 ) - v£ 



1 



u(u - e k T 3 ) - v 2 . 
Evaluating the commutator 

[V+VV^a] = 2(ui k Vv 2 k - «2 k V«ik)n, 
we may then write 
4e 2 T 



(156) 



(157) 



(158) 



(159) 



^(«ikV« 2k - « 2k V Vlk ) 2 Tr[G( K ) 2 ] , (160) 



where for convenience, we have assumed an isotropic stiff- 
ness Q^b = Q 5 a b- Taking the zero-temperature limit, 
replacing iu> n — ► iuj and 



to' 



(161) 



2 c 
e P s = 



^(fikVu 2 k - W2kVwik) z 

k 

5 du> 1 



(162) 



Carrying out the final frequency integral yields the fol- 
lowing expression for the composite stiffness 

= e*y(vjv 1 \ 2 (* lk V* 2k - ^kV^k) 2 
d^y Vk+ J [e k + (2 Uk+ ) 2 ]^ 

This result is quite fascinating, because it confirms that 
the composite pair condensate has a phase stiffness, even 
in the absence of an underlying Fermi surface. In a con- 
ventional superconductor, the scale of the stiffness is de- 
termined by the Fermi energy e 2 p s ~ e 2 ep/o, 2 . Here the 
size of the stiffness 



2 c 
e Ps 



e 2 N(0) 




e 2 T c /a 2 



(164) 



is determined by the condensation energy. This is a clas- 
sic example of "local pair" condensation. 

Away from half filling, the superfluid stiffness contains 
contributions associated with both the paired conduction 
electrons, and the composite pairs. To make a clean di- 
vision of the stiffness into two components, one needs 
to worry about the various cross-terms that appear in 
the stiffness such as ( |153| ) , which do not obviously vanish 
away from half filling. However, a key observation is that 
these terms are odd functions of A, so they are guaran- 
teed guaranteed to vanish provided that the physics is 
particle-hole symmetric about A = 0. With this proviso, 
we can divide the superfluid stiffness into two terms 



Q = Q 



BCS 



Q c 



(165) 



where 



T ■ 



[iG(K),T 3 ] 



then yields 



Qab = |Tr[[G c/ V a V,r 3 ][G c/ V fc V,r 3 ] +H.c.]. (166) 

are the BCS and composite pair contributions to the stiff- 
ness. 



XI. DISCUSSION 

In this discussion, we should like to address the results 
of this paper on two fronts: 

• Alternative theoretical approaches to test and con- 
firm the presence of a co-operative Kondo effect. 

• Applications to experiments and the theory of real 
heavy fermion systems. 
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A. Alternative Theoretical Approaches 

The main effort of this paper has been to establish 
the link between co-operative channel interference in the 
Kondo lattice and composite pairing. Although our key 
result- the relation between the gauge invariant interfer- 
ence term and composite order 



V f 2 Vi 



J1J2 



Ft A 

-At Ft 



(167) 



does not depend on approximations, the notion that this 
interference term can develop long-range order relies on 
various mean-field approximations. There are however 
reasons to be confident in the mean-field solution. First, 
there are no obvious competing instabilities, such as mag- 
netism. By turning on a J2, one does not drive the system 
closer to the antiferromagnetic instability, but rather, 
simply activates another source of Kondo screening. This 
should be contrasted with the situation in the alternative 
model of spin-fluctuation mediated pairing. Here, to at- 
tain a transition temperature that is comparable with the 
heavy fermion band-width places the model close to an 
antiferromagnetic instability, where the competing effects 
of magnetism make the mean-field theory potentially un- 
reliable. 

Nevertheless, the use of mean-field methods inevitably 
raises questions about our work which motivates us to 
seek alternative methods to verify the key results. It 
may be possible to precisely verify our results in both 
finite size calculations, and in the exactly solvable limit 
of infinite dimensions. Finite size studies on our model 
may be facilitated by treating the model as a ladder com- 
pound and by using the strong-coupling limit so as to 
completely eliminate the possibility of antiferromagnetic 
instabilities. 

To formulate our model in a form that is tractable to 
an exact infinite dimensional study, rather than using ex- 
plicit form factors in the Kondo interaction, it is better 
to start with a channel conserving two-channel Kondo 
lattice, to which a term which destroys channel conser- 
vation is then added to activate the channel interference. 
Suppose one starts out with a two-channel Kondo model, 
with perfect channel conservation: 

h 



H c = 22 ekctrkfjcrkfj + ^ Jrc^rj<rcrj ■ Sj . (168) 



kfer 



where 



J rja = N s 2^ c t rkCT£ 



-ik-R, 



(r = 1, 2) and 

e = —It ^ (cos ki) — jj,. 



(169) 



(170) 



This model has a perfect U(l) channel symmetry. Sup- 
pose one now adds a term to the hopping H — > H c + H', 
where 



H' = -2A 



;=i,< 



(-l) i COsfc ( (c t kl(T Ck2 ( T +ct k2CT C k la), (171) 



is a hopping term with "d-wave" symmetry that mixes 
the different channels, breaking the U(l) channel conser- 
vation symmetry. This contains no no non-local interac- 
tions, and is thus ideal for a large d-treatmcnt. It is easy 
to see that in the limit of large A, it is equivalent to a 
one band model, with two orthogonal form factors 



$ik 

$2k 



1, 



(172) 



To see this, note that when A is non-zero, the band splits 
into two components with energies 



F k ± = e k ± 2A 



(-l)'cosfcz 



:=i,3, 



(173) 



where 



kcr 

1 



Oko- = — ^(Ckla + < f ) 2kCk2 CT ), 

1 

bka = ~^( Ckla ~ <f> 2kCk2 C r)- 



(174) 



At large A, one can project out all terms involving the 
upper band by re-writing the Hamiltonian in terms of 
the a and b creation operators, then all terms involving 
the upper band creation or annihilation operators ak CT or 
a^kcr- When we do this, the interaction becomes 



Hi = tt^-Y] Jk.k'^Vfek • V'( k - k ')-*; 



k.k' 



Jk,k' = (Jl*lk*lk' + ^2*2k*2k')> 



(175) 



corresponding to a two-channel Kondo model with two 
orthogonal matrix elements. 

The mean-field theory for this model predicts that 
that as soon as A becomes finite, channel interference 
will drive the system into a composite paired state, even 
when Ji >> J 2 . This phenomenon should extend all the 
way out to infinite dimensions, where a precise dynam- 
ical mean-field theory treatment of the model becomes 
possible. 



B. Possible consequences of co-operative Kondo 
behavior for Heavy Fermion compounds 

We conclude with a brief discussion of the general con- 
sequences of co-operative channel interference in heavy 
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fermion systems. Our paper has focussed on the super- 
conducting aspects of this problem. Here we should like 
to put the problem in a more general perspective. 

Assuming it becomes possible to verify the theoretical 
soundness of the co-operative Kondo effect, how could 
the theory be tested experimentally? 

• Our model suggests a rather intimate relation be- 
tween the local quantum chemistry of the heavy 
fermion ion, and the gap symmetry of the order 
parameter. In heavy fermion compounds, one of 
the scattering channels is an f-channel. Since the 
two channels must have the same parity, the second 
channel is in all likelihood another f-channel, or a 
p-channel: 

— * / ® / ■ Candidates: Non-Kramers ion, 
e.g URu2Si 2 , UBei 3 . 

— o / <8> p ■ Candidates: UBe\s, UPt 3 and 
Cerium systems, close to quantum critical 
point. 

The first possibility will occur if the Kondo effect 
involves a non-Kramer's magnetic ion. For exam- 
ple, in the case of URU2S12, there is strong cir- 
cumstantial evidence that the single-ion physics is 
dominated by a Kondo effect with a non-Kramers 
magnetic doublet, The form factors for the two f- 
channels in a tetragonal crystal field are known, 
and place strong constraints on the symmetry of 
the putative composite order. In the second case, 
the number of available p-channels is small and the 
local quantum chemistry will determine the most 
likely channel for the co-operative pairing process. 
For example, in hexagonal UPt 3 , the most likely 
second-channel is the p z orbital, which would ex- 
plain the presence of the node in the basal plane. In 
principle, cubic U Be\^ could belong to either cate- 
gory, as this system may also have a non-Kramer's 
ground-state. If however, the driving force derives 
from a p-channel, molecular orbital theory dictates 
that the most likely second-channel is a p-wave 
state with normal orientated along the cube diag- 
onals, such as the 111 direction. A gap node nor- 
mal to this direction could be detected-using careful 
transverse ultrasound measurements. E_l 

• One of the paradoxical features of heavy fermion 
superconductors, is that their large entropy of con- 
densation suggests a large superconducting order 
parameter. Yet Josephson tunneling with a con- 
ventional superconductor has not, to date been 
achieved. Although composite and normal pairs 
co-exist side-by-side in our hypothetical super- 
conductor, the predominantly composite character 
of the order parameter may help explain why it 
has proven so interminably difficult to carry out 
Josephson tunneling into these systems. One way 



to enhance the Josephson current may be to intro- 
duce rare earth, or actinide spins into the tunnel 
junction. Josephson tunneling between a conven- 
tional, and composite paired superconductor re- 
quires that the addition of a pair is accompanied 
by a spin-flip. Spin fluctuations of the local mo- 
ments in the junction may help to catalyse this co- 
operative process. This is a possibility currently 
under investigation. 

We should like to end with a short note about the non- 
superconducting aspects of the composite Kondo effect. 
In our key identity 



V^Vi 



J1J2 



F A 

-At Ft 



(176) 



we have the possibility of finite diagonal components 
F^O due to co-operative interference. Unlike composite 
pairing, such instabilities will require J\ and J2 to be of 
comparable size. There are, to our knowledge two good 
candidates for this kind of phenomenon: 

• Orbital magnetism in U Ru2S%2- As mentioned 
above, this material is a naturally occurring two- 
channel Kondo lattice, but with strong spin-orbit 
coupling. One of the long-standing mysteries of 
this compound, is the appearance of an uniden- 
tified magnetic state at 17K, with a large order 
parameter which appears to break time-reversal 
symmetrii_iu.it without producing a large magnetic 
moment P 1 " 32 ! One possible way to account for this, 
is to suppose that the two channels in this com- 
pound give rise to a complex order parameter 



F(x) = F e iQ ' x , 



(177) 



where Q is commensurate with the lattice. Just 
as superconducting composite order coexists with a 
weak BCS order parameter, orbital composite order 
will co-exist with a weak orbital moment. Spin- 
orbit coupling will then generate a weak magnetic 
moment. 

• Ultra- narrow gap Kondo insulators, CeRhSb and 
CeNiSn. These Kondo insulators appear to de- 
velop gap nodes iu their tiny hybridization gap. 
In a recent paper J_3 we have pointed out that this 
kind of behavior would arise from the suppression 
of shape fluctuations, which gives rise to three or- 
bital scattering channels in which the Kondo effect 
can take place. The resulting interference between 
the three orbital channels is found to spontaneously 
generate a crystal field environment that gives rise 
to a Kondo "insulator" with gap nodes. 

These are both areas of active investigation, which lie 
outside the scope of this paper. 
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APPENDIX A. 



In this section we apply the SU(2) decoupling scheme 
originally developed by Affleck et al. for the Heisenberg 
model, to the two-channel Kondo problem. This deriva- 
tion is closer in spirit to the original work by Affleck et 
al, and differs in detail from the later work by Andrei and 
Coleman. The approach of Affleck et al is more explic- 
itly SU(2) symmetric and can be naturally extended to 
include source terms that couple to the composite order 
parameter. There are two distinct differences between 
the approaches: 

• The integration measure over the SU(2) field is 
"fiat" 

• The Gaussian coefficient of Tr[VVVr] is now 1/J 
rather than 1/2 J, as it was in the earlier work by 
Andrei and Coleman. 

The difference in measure leads to differences in the 
fluctuations around the mean-field theory, and the mean- 
field expressions for the Kondo temperature obtained in 
the two methods actually differ by a factor of two in 
the exponential. We have chosen the approach of Affleck 
et al because it gives us a much cleaner an symmetric 
derivation of the final results. 

The objective of this section is to show how the inter- 
action between a localized moment S and the electron 
local spin density or = i/Atrip can be decoupled in terms 
of a fluctuating SU(2) field 

Jr(S ■ a r - ±) -> [/ f V r # + H.c] + Tl[V f Vr] (its) 

■Jr 

where Vr is directly proportional to an SU(2) matrix gr 

r 

Vr - iV^gr 



V A 

A* -V* 



(179) 



For clarity, all site indices j are omitted from this deriva- 
tion, but are readily restored later. 

Following earlier work, we introduce the following ma- 
trix fermions 



/T h 

f\ -A 



4i 



-4 T 



(180) 



By taking the product of these matrix operators with 
their Hermitian conjugates, we find that 



*r*r = 1 + £ T • ffr 
P T = 1 + 2<r T • S 



(181) 



where <Xr = ip^r&ipr is the electron spin density and ct t 
denotes the transpose of the Pauli spin matrix. Transfor- 
mations acting to the right of T, T — > Th correspond to 
physical rotations of the local moment. Transformations 
acting to the left of T , T — > gT correspond to the local 
SU(2) transformation, under which the spin operator is 
explicitly invariant: 



(182) 



Multiplying the two equations (181) together, and taking 
the trace we obtain 



(183) 



Anti commuting the conduction electron operator to 
the left through the trace, we then find that 



Jr(<rr-S-i 



Jr 



T±[UpJ v \. 



where 



(184) 



(185) 



is an anti-unitary matrix and 



(186) 



Notice that if we expand the above interaction, we obtain 

Hi = -'^[a\a v + a r a r + & r &r + b r bl] (187) 

showing that it has been decoupled simultaneously in the 
particle-hole, and Cooper channels. 

We now apply a Hubbard-Stratonovich procedure to 
this expression. Formally, we first convert each of the 
fcrmionic operators in the interaction to Grassman vari- 
ables inside a path integral. On each time slice we write 



-AtHi 



D[Vr,V^rV 



-AtHi [Vr.V+r] 



where we have transformed 



Hi = ~Ti[U^U r ] -> F/[Vr,V + r] 



(188) 



(189) 
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and 

Ulc^Ur] + Tr[E/ r V r ]} + ^-Tr[V r V r ]. (190) 

A priori, V is a two by two complex matrix. However, 
if we divide it up into the sum of a unitary and an anti- 
unitary matrix, we find that only the former completely 
decouples. The residual part of V is completely anti- 
unitary, and has the form 



V r = 



V 

A* 



A 

-V* 



(191) 



where there are only two independent complex parame- 
ters. This is a significant simplification. Notice that Vr 
is directly proportional to an SU(2) matrix: 



V r = iV r g r , (V a r = vW+W) ■ 



(192) 



The measure of integration for each time-slice is then 
simply 



D[V r , V f r ] = dV r dV^ dA r eiA r . 



As our final step, we now reduce the decoupled inter- 
action to a more manageable two-component notation. 
Writing 



fir = (^r-[,ipri), 
then Hj reduces to the form, 



H, 



[/tV r # + he] + — Tr[V r V r ] 
Jr 



(194) 



(195) 



which is the form quoted in the main text. 



APPENDIX B. 



The purpose of this section is to establish the direct 
relationship 



V f 2 Vi = - 



J1J2 



Ft A 

-At Ft 



where 



F = tfainfa ■ S, 



(196) 



(197) 



represent the single composite order in the particle-hole, 
and particle-particle channels respectively. 

In order to establish this identity, we introduce a source 
term into the Lagrangian which couples to the gauge in- 
variant matrix product V2V1, writing 



Hi = ^|[/ f V r # + he] + i-Tr[V r V r ]| 



TrD^Via + h.c] 



(198) 



where the source term a = a + ia-T is a unitary matrix, 
with four real coefficients. 

We shall now invert the Hubbard Stratonovich trans- 
formation, with the source terms in place. We begin by 
rewriting the Gaussian term in the interaction to obtain 

H i = Y, iP^r^r] + H + Tr^Vr'mr/r] (199) 



where 



TOFT' 



77 a 



(200) 



When we carry out the Gaussian integral over Vr, the 
transformed Hamiltonian now becomes 



H r 



1 



Tr[[/t r £/ r ,( m -Vr], 



(201) 



(193) where 



J1J2 



(i-H 2 JiJ 2 ) 

When we expand this to linear order in a we obtain 

JiJ: 

r 



(202) 



Hi = - ^Tr[C/t r[/r ] + ^l T r[U\U ia + he] (203) 



Inserting T^T = 1 + 2<x T • S into this expression, and 
making the observation that 

Tr [0^2*1 + he] = (204) 

we can rewrite 

Tr[C/t 2 [/ ia + he] = 2Tr[a* 2 (o- T • S)* t 1 + he] (205) 

so that the final form of the interaction, with the source 
term is 



ff ; = -^iTr[[/t r t/ r ] 



-Tr[a* 2 (<r T • S^x + hc] 



(206) 



comparing coefficients of a in ( |198| ) with ( |206| ), we obtain 
the following identity 



VfoVx 



J1J2 



(207) 



This expresses, in a compact form, the relationship be- 
tween the inter-channel interference and the composite 
order. To complete the job, we now expand the right- 
hand side. We first write 
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[<& 2 (<7 T • S)&i] = -[*2(<T y V* y ■ 8)*\] 



(208) 



XAA 



where we have replaced <x T = — (a y <7a y ). To make the 
expansion, it it is convenient to write 





"o 


" 







" 




\G K 





T 3. 


Qk 





T 3_ 


) 



(215) 



At half-filling, the electron propagator can be written 



^2 = 



(209) 



where is the row-spinor formed by taking the trans- 
pose of the column spinor tf>2, and ip* — (i/'^i) T is the 
column spinor formed by taking the transpose of ip 'i- 
Multiplying out the matrices, we obtain 



JiJ: 



■Vf 2 Vi 



(jycrayipl ip2 ' (-io-ycr)ipi 

2 {ci<J y )i>t i>2^<ri>i 

ift 2 (<Ti(T v )lpt ^(Tlpx 

F A 

-At Ft 



Q{k) = 

where t>k = \/ v( k + v 
G(k) 

and 



ujG —iv^Gg^ 
ivkgG g(uj - ekT 3 )Gg J< 



T- 

2k- 



lu(lj - e k T 3 ) - vl 



■ S 



g = [Vlk - lV2kT 2 \ 



(210) 



The Green function Q has poles at iE^ where 



Substituting (210) into (207)) we obtain the quoted re- 
sult. 
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The purpose of this section is to evaluate the suscep- 



(v = ±) 



(216) 



(217) 



(218) 



(219) 



Expanding the electron propagator about its poles, we 
write 



73 



r,=± 



G v (k) 



u - E kv r 3 



(220) 



tibilities associated with the expansion of the mean-field i nser ting this into the full propagator, we can write it in 
Free energy about the pure composite paired state ^e form 



F = F °~2 [Xm*) 2 + 2 ^a<^A + xxx(SX) 2 ] . (211) 

By integrating over the Gaussian fluctuations in A to im- 
pose the constraint on the f-charge, we can use these 
susceptibilities to compute the physical charge suscepti- 
bility 



ci v G v -ick n s kv G n T 3 g^ 
ic^s^gTzG^ s^ v gG v g^ 



m = J2 



where 



Xc = 



Xaa 



(212) 



Ck); — 



isk v gkT 3 



(221) 



(222) 



To compute the susceptibilities we expand the Hamil- 
tonian about the half-filled state, 



is the eigenvector corresponding to the quasiparticle with 
energy E^. The quantities 



H = H — Sfi 



t 3 




+ A 





T 3 



The electronic part of the Free energy is given by 
F e = -T^Trln[io; n - H(k)} 

K 

Expanding this to second-order then gives 



(213) 



(214) 



■ ki/ 



-ki; 



1 



E^ — ^k-j) 2 

2 E^ — £k 1 

k,? E^ — i?k- t; 2 

- _ 



1 + 77 

1-7? 



fk 



£k 



(223) 
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describe the admixture between the conduction, and f- 
electrons. The matrix elements of the charge operators 
appearing inside the susceptibilities are 





T 3 



Ck, f 



T 3 
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o o 

o T 3 



f 2 

= s^Skri' (cos 0k r 3 + sin T 2 ) (224) 



where Ck = «k_/ u k an d = 2vikU2kA>k- The expres- 
sions for the susceptibilities can now be written 



Xw = ~ T E Tr 

re, 77,77' 

>^a = T E Tr 

K, 77,77' L 

-T ]T Tr 



XAA 



G^TzGrjiTzC^ + G^GrfTiSl- 



Xb 



E 



/ 2t»lk^2k \ ' 



[(e k )2 + 4 ^]3/2 ^ „2 J 



(5+4). (226) 



Since xaa > Xmm> and IX/*a| < Xw tne tne nnal result for 
the charge susceptibility is then guaranteed to be positive 
when V1V2 =/= 0: 



Xc = Xmm 



(Xma) 2 
Xm + Xb 



> 0, (wit* # 0). (227) 



where we denote G v = G v (k) and vanishing cross-terms 
between T3 and t 2 have been dropped. We now evaluate 
the Matsubara sums in these expressions, and take the 
zero-temperature limit. Key results that we use are 



T £ Tr 

UJ n ,Tj' 



GnTaGrfTa 



2 2 _ oV" -/W f kr l -2-2 

-T7V ^Z^ p, _ p, ,0' 



2c 2 s 2 



where fu v = l/(e' 3 • Ek ' , + 1) denotes the Fermi function. 
Similarly, 



T E Tr 



GnTzG^Tz 



(cs) v (cs) v , 



2c 2 s 2 



T E Tr 



„2 „2 



Lastly, there is one anomalous term, 

G I) T 2 G r) /T2 



2c 2 s 2 



T Tr 



2 2 
s r) S r]' 



Putting these results together, we obtain 



Xma» 
Xa*a 



E 

k 

E 



4v 2 



k ( £ k + 4 «k) 3 



4vt 



k + 



XAA = J2 



k (4 + K) f 



2 1 e k / 2vi k v 2 k \ 2 



k ^v / 4 + KV vl J 

We can also rewrite xaa in the following form 

XAA = Xw + Xb 



(225) 
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